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The problem of entanglement produced by an arbitrary operator is formulated and a related mea- 
sure of entanglement production is introduced. This measure of entanglement production satisfies 
all properties natural for such a characteristic. A particular case is the entanglement produced by 
a density operator or a density matrix. The suggested measure is valid for operations over pure 
states as well as over mixed states, for equilibrium as well as nonequilibrium processes. Systems of 
arbitrary nature can be treated, described either by field operators, spin operators, or any other kind 
of operators, which is realized by constructing generalized density matrices. The interplay between 
' entanglement production and phase transitions in statistical systems is analysed by the examples 

of Bose-Einstein condensation, superconducting transition, and magnetic transitions. The relation 
between the measure of entanglement production and order indices is analysed. 
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I. INTRODUCTION 



> 

in 

I Entanglement is the term used by Schrodinger with respect to the superposition principle applied to composite 
. systems [1-3]. If two quantum particles have interacted, their state cannot be presented as a tensor product of single- 
particle states, but it is entangled, being a superposition of such products. The notion of entanglement is now at 
. the heart of such interrelated intriguing problems as quantum measurement, quantum information processing, and 
■~»>^ ' quantum computing, which have been expounded in several books and reviews [4-12]. 

, In the literature, one distinguishes between the entanglement of quantum states and the entanglement produced by 
^Y^' quantum operators over unentangled wave functions. Both these types of entanglement, to be well defined, require 
[ the knowledge of quantifying characteristics. In the present paper, as follows from its title, we shall consider only 
^ ■ the second type of entanglement produced by quantum operations. In quantifying the latter, one usually envokes the 
^ ' characteristics measuring the entanglement of states. Because of this, it is useful to briefly mention the problem of 
quantifying the entanglement of quantum states, which, at the same time, would make clearer the difference between 
' these two types of entanglement. 
. ^ , Quantifying entanglement of quantum states, one usually deals with bipartite systems. Several ways of measuring 
• entanglement in such systems have been suggested, the methods being based on the notions of either reduced or 
' relative entropy. Thus, mutual information is a linear combination of the von Neumann reduced entropies [13,14]. 
- - 1 The reduced entropies themselves define the entropy of entanglement, which serves as a measure for entanglement of 
formation [15,16]. Another measure of entanglement is defined by minimizing over disentangled states the Kullback- 
Leibler distance, yielding the relative entropy of entanglement [17-19]. A measure, not envolving the notion of entropy, 
could be introduced as the number of maximally entangled pairs that can be purified from a given state, which results 
in entanglement of distillation [15]. However, this measure depends on the particular process of purification, and it is 
not clear how to compute it in an efhcient and unique way [17-19]. 

Entanglement of formation was also employed in the attempt of characterizing the entanglement for mixed states 
[20], but it was concluded that this way did not uniquely define mixed-state entanglement [20,21]. A suggestion for 
measuring covariance entanglement in bipartite systems by means of correlation functions squared was advanced in 
[22]. However, such a measure, as is accepted by the authors themselves, does not possess all necessary properties 
[8-12] to be considered as really a measure. One often connects the existence of entanglement with the violation of 
the Bell inequalities, which can be formulated for bipartite systems with both orthogonal and nonorthogonal states 
[23]. This point of view is based on the Gisin theorem [24], according to which any pure entangled state of two 
particles violates a Bell inequality for two-particle correlation functions. Nevertheless, there exist pure entangled 
N > 2 qubit states that do not violate any Bell inequality for TV-particle correlation functions [25]. It seems that 
correlation functions do not provide the best tool for characterizing entanglement. At present, there is no such a 
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general definition of entanglement measure that would be valid for bipartite as well as multipartite systems, for pure 
as well as mixed states, for cqiiilibrium as well as nonequlilibrium processes. 

Another problem is to quantify the entanglement produced by quantum operations on a given set of disentangled 
functions. The so produced entanglement is termed as entangling power, entanglement capacity, entanglement of 
evolution, entanglement generation, or entanglement production [26-30]. In what follows, we shall employ the term 
entanglement production as the most closely related to the meaning of this notion, being the entanglement produced 
by an operator. To quantify this type of entanglement, one usually resorts to the combination of measures defined 
for the entanglement of states. Since the latter are well defined only for pure bipartivc systems, the entanglement 
production is also usually considered for such systems. 

The aim of the present paper is to introduce a general measure of entanglement production, which could quantify the 
amount of entanglement produced by an arbitrary operator on a given disentangled set. The suggested approach can 
be applied to operators of any nature and to any physical systems, whether pure or mixed, bipartite or multipartite, 
equilibrium or not. Being justified for arbitrary operators, the approach can straightforwardly be applied to a 
particular kind of operators, as statistical or density operators. 



II. DISENTANGLED AND ENTANGLED FUNCTIONS 



First of all, it is necessary to give correct mathematical definitions for the notions of entangled or, conversely, 
disentangled functions and to specify the notations to be employed in what follows. From the very beginning, a 
multipartite composite system is kept in mind, consisting of an arbitrary number of subsystems enumerated by the 
index i = 1,2, ... ,p, where p = 1,2,... can be any integer. Subsystems are treated as indivisible parts, because of 
which they could equivalently be called particles. These can be distinguishable or indistinguishable, which does not 
change the general mathematical structure, provided an appropriate labelling of partite states is employed, e.g., by 
means of collective mode labels or occupation numbers [31-33]. 

The space of single-partite states for each i-part is presented by the Hilbert space 

Hi = Z{\ni >} , (1) 

being a closed linear envelope of a basis {\ni >} composed of orthonormalized single-partite states [rij >. Hence, any 
vector ipi e Hi can be expanded over the basis {[n, >} as 

n. 

The nature of the labels rij here is of no importance. The vector norm 

in Hi is defined through the associated scalar product {ipi, ipi). 

It is worth emphasizing that |n > does not compulsorily mean a quantum state of a given physical particle. This 
case is not excluded if particles are distinguishable. But, if one deals with a system of indistinguishable particles, \n > 
should be understood as a single-particle mode. 

The space of composite-system states is given by the p-fold tensor product 

H = > (3) 

which is identified [34] with the closed linear envelope 

H = C{\ni...np>} (4) 
over a normalized p-particle basis {jni ... rip >}. The latter may be written as the tensor product 

jni .. .rip >= 0^^Jnj > (5) 
of the single-particle basis states. Any function ip gH can be presented as the sum 

'P=^ Cni...n^\ni...np> (6) 
{ni} 
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over the multiparticle basis {\ni . . .Up >}, where {rij} = ni,n2, . . . ,np. The vector norm ||<^||-h = ■\/(Jp7v) '^^ is 
generated by the related scalar product {ip, (p). 

Note that it is not compulsory to deal with the total tensor- product space (3). In some cases, because of physical 
restrictions, additional selection rules may be superimposed on the admissible states of Ti. This is, e.g., the case 
of systems composed of identical particles, whose quantum states are to be either symmetrized or antisymmetrized, 
according to whether the particles are bosons or fermions. Then the space of admissible states is reduced to a subspace 
of space (3). In what follows, we shall keep in mind the possibility of such additional restriction rules. For short, 
we shall continue denoting by 7i a subspace of the space structure (3), implying that all necessary selection rules are 
taken into account. In the presence of the latter, the structure (3) is termed the incomplete tensor product [35 39]. 

The multiparticle space (3) contains the whole variety of different states, from which one should separate out 
disentangled states or factor states that are presentable as the tensor products (8)f=i<^i. A collection {/} of all possible 
factor states / forms the disentangled set 

V = {f = , ifi e Hi} , (7) 

which is a subset of the Hilbert space (3). A function f gV cH can be written as 

where expansion (2) is taken into account. For any two vectors / and /' from V, the scalar product is 

{fJ') = fli^i,<p'i)- (9) 

i=l 

This generates in T> the vector norm 

ii/ih=nii^.ik- (10) 

i=l 

The compliment T-L\'D to V forms the set of entangled states. 

To illustrate in an explicit way the principal difference of disentangled states from entangled ones, let us consider 
a bipartite system, with p = 2, and let the single-particle states be two-dimensional. Then, writing, for compactness, 
1 and 2 instead of ni and 712, for state (6), we have 

(/3 = Cii|ll> -t-Ci2|12> -FC2l|21> -hC22|22>, 

while for the factor state (8), we get 

/ = ai6i|ll> +ai62|12> + a26i|21 > + 02^2122 > . 

As is evident, the state ip £ Ji is more general than / e V. The space Ji contains the entangled states, such as 
C12II2 > -|-C2i|21 > or Cii|ll > -I-C22I22 > that in no way can be reduced to the factor states / € P. In general, no 
entangled state can be presented as a product of single-particle states. 



III. MEASURE OF ENTANGLEMENT PRODUCTION 



The multiparticle space (3) contains both entangled and factor states. For the time being, we follow the abstract 
terminology of Sec. II, implying no physical applications that will be treated later. An abstract mathematical level 
of consideration provides the best way for making transparent what actually is entanglement production and how to 

measure it. 

The term entanglement production as such means an action that transforms disentangled states into entangled ones. 
A transformation can be described by the action of an operator. Thus, one may investigate entanglement produced by 
different operators. Let A be an arbitrary linear bounded operator acting on the tensor-product space W. A complete 
theory of linear bounded operators, defined on the tensor-product spaces, can be found in Refs. [35-39]. 

The norm of a linear operator A onTi. can be given by 

||A||„= sup \\A^\\h= sup \{^,Aip')\. (11) 

\\v\\n=l \M\-H=\\v'\\n=l 
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For a bounded operator, the norm is finite. If the operator A is Hermitian, then 

||A||„= sup \{^,A^)\ {A+=A). 

\\v\\n=i 

Let us introduce the projector Vd, which projects the total space (3) onto its subset given by the disentangled set 
(7), so that 

VvH = V , (12) 
with the standard properties of projecting operators 

Vl = Vv, V+ = Vv, \\Vv\\n = l- 

The projector in Eq. (12) is nonlinear, therefore the equality for its norm has to be understood as a definition. And 
let us define the norm of A on 2? as 

\\A\\'d = \\VvAVv\\h ■ (13) 
This, in view of the structure of the disentangled set (7), can be presented as 

||A||^= sup \{f,Af)\. (14) 

||/||i.=||/'||i.=i 

An operator A, acting on / G I?, generally, transforms / to an G 7i. That is, an operator A^ in general, entangles 
the factor states. The operator A, having the property AV C H \ I? will be termed entangling operator. 

Similarly to the existence of entangled and factor states, there exist entangling operators and nonentangling ones. 
The latter should, clearly, have the structure of a direct product '^^^iA\ of single-particle operators A\ acting on 
TLi. Lot the algebra of all linear bomidcd operators on 7i be denoted hy A= {A}. And let us separate out from this 
algebra a subset = {A'^} C ^ of nonentangling, or product, operators having the structure of a product (Si^^iA\. 
Thus, by construction, 

AVcn , A®X> c V . (15) 

Analogously to the projection of H onto V by means of the projector Vt> given in Eq. (12), we may denote the 
reduction of A to A'^ with the help of a projector V^, such that 

A^=r^{A), (16) 

which is a superoperator acting on the Hilbert space of linear operators, with the properties 

^|=^», H = 'P^^ ll^»lb = i 

being valid. The equality for the norm has to be understood as a definition. An explicit construction of a product 
operator A® c A^, associated with a given operator A c A, can be done as follows. Let us define a reduced 
single-particle operator A\ on Hi as 

A{ = Ci TV{„,^,} A , (17) 

where the trace runs over all Tij except the case j = i. The set of constants Cj is chosen so that to satisfy the 
normalization condition 

Tr^^ A = Trx. A® . (18) 

In this way, we obtain the product operator 

T^v <=i A\ 

where 



Ci 4 . (19) 



Tri, 4 = i]Tr„, A[ . 



i=l 
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Thus, from a given multiparticle space (3), it is possible to separate out the disentangled set (7) and, similarly, for 
an arbitrary operator A, one can put into correspondence the product operator (19). That is, there exist entangled 
and disentangled states as well as entangling and nonentangling operators. 

Let us be interested in the entanglement produced by an operator A onTi. How could we measure the resulting 
entanglement production? It would seem natural that we should somehow compare the actions of the given operator A 
and its nonentangling counterpart A® . But what quantity should be defined for this purpose? Following the common 
ideology, one should construct a sort of entropy, comparing, say, Tr-^AlnA and Ti-pA'^ In A®. However, this is not 
the best way. The pivotal idea, we shall follow, is the observation that the norm of an operator characterizes a kind 
of order associated with this operator [40]. For instance, invoking the norm and trace of reduced density operators, 
or reduced density matrices, it is possible to define the density order indices [41]. Generalizing the latter, one can 
define the operator order indices [40] for arbitrary operators. These order indices provide a complete classification for 
different types of order, long-range and mid-range, off-diagonal and diagonal, because of which they are applicable 
for describing both phase transitions and crossover phenomena. 

Entanglement, in some sense, is also a characteristic of order (or disorder). Hence it should be well characterized 
by an operator norm. To this end, we introduce the measure of entanglement production 

generated by the given operator A. Here the logarithm can be taken with respect to any base that would be convenient, 
e.g., to the base 2. 

As follows from definition (20), the entanglement-production measure quantifies the amount of entanglement pro- 
duced by an operator A over a set V, because of which this measure should, in general, be denoted as e{A, V). When 
working with a fixed set T), one may, for simplicity, shorten the notation writing e{A) = e{A,'D). Quantity (20) 
satisfies all natural properties that are compulsory for being really a measure: 

1. Measure is semipositive. 

This is evident by construction, since 

ll^^lh = 11^0(^)112. <ll^lh, 

because of which 

s{A) > . (21) 

2. Measure is continuous. 

This implies the following. Assume that for any operator A of H there exists a family {A{t)} of operators A{t) 
parameterized with f e R, so that A{t) ^ A as t ^ in the sense of the norm convergence 

\\A{t)\\T,^\\A\\T, (t-0). 

If so, then for the measure (20) one has 

e{A{t)) ^ e{A) {t ^ 0) . (22) 

The nature of the parameter t can be arbitrary. In physical applications, this may be time, temperature, density, 
interaction parameters, and so on. 

3. Measure is zero for a nonentangling operator. 

A nonentangling operator has the form of A®. As is obvious, for A = A®, 



= . (23) 

In particular, there is no self-entanglement of a single-particle system, when A = Ai= A®. 
4. Measure is additive. 

Let now {A^} be a set of copies of A, with v = 1,2, . . ., so that A = ®vAv. Such a case may happen, e.g., in 
treating heterophase systems [42]. Then A® = ®uA® and one gets 
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Prom here, it follows that 

e A,) = e{A,) ■ (24) 
5. Measure is invariant under local unitary operations. 

Such operations are described by a set {Ui} of unitary operators Ui on Hi, with U^Ui = 1 and i = 1,2, . . .,p. 
Using the properties of the operator norm and trace, it is easy to show that 

ei^UU^A^U Ui)=eiA). (25) 

Thus, the entanglement-production measure (20) can be introduced for any operator. The requirement that the 
latter be bounded can be relaxed in the following way. Assume that A defined on H is unbounded. Introduce a 
restricted space TZn such that the norm ||A||7^„ is finite and H can be treated as the inductive limit TZn ^ H 
as N ^ oo. Here N is not necessarily the number of particles, but may be any labelling number. Following the 
procedure described above, one can separate out from the space TZn the disentangled set Pat, whose inductive limit 
isV = limjv->oo f jv- Then the entanglement measure, instead of Eq. (20), can be defined as 

e{A)^ lim log 



In particular, this can correspond to the thermodynamic limit. 

IV. PURE STATE ENTANGLEMENT 

Entanglement production, as is explained above, can be generated by any operator. In physical applications, one 
often implies that this is due to the action of a von Neumann operator. There are two types of density operators, 
pure-state and mixed-state. A density operator ppf for an ensemble of N particles is termed a pure-state operator 
when it is an idcmpotcnt operator, such that = p^j. Below, we shall consider several examples of calculating the 
entanglement-production measure (20) generated by pure-state density operators. 

A. Einstein-Podolsky-Rosen states 

These states provide a classical example of entangled states. They correspond to a bipartite system composed of 
two-dimensional parts, so that such a state can be presented as 

|£;Pi?> = ^ (|12 > ± |21 >) . (26) 
V2 

The related density operator is 

PEPR = \EPR >< EPR\ {TrnPEPR = 1) , (27) 
with W = Wi (g) Reduced single-particle density operators, according to definition (17), have the form 

p\= i(|l><l| + |2><2|), 

being given on Hi, with i = 1,2. The product operator (19) is pf = /3} ig) p^. Calculating the norm (13), we have 
for the density operator (27) 



while for the product operator. 



\Pepr\\-d = llPillw. = 2 



WPeprIw = WPiW-Hi = ■ 
Thence, the entanglement-production measure (20) is 

e{pEPR) =log2. (28) 
Assuming the logarithm to the base 2, one has £{pepr) = 1- 
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B. Bell states 



Similarly to the previous case, it is easy to find the entanglement-production measure for the density operator 
formed by the Bell states 

|5>=^(|11> ±|22>), (29) 

so that the corresponding density operator is 

PB = \B >< B\ (Tr« PB = 1) . (30) 
For the entanglement-production measure (20), this gives 

e{pB) = \og2. (31) 

Keeping in mind the logarithm to the base 2 yields £(/0b) = 1. 

Note that the Bell states can be treated as a two-particle generalization of the single-particle Schrodinger cat state. 



C. Greenberger-Horne-Zeilinger states 

Such states, having the form 

\GHZ > = (|11 . . . 1 > ± |22 . . . 2 >) , (32) 

can be considered as an iV-particle generalization of the single-particle Schrodinger cat state [43,44]. The associated 
density operator is 

PGHZ = \GHZ >< GHZ\ {Tr-HPGHZ = 1) • (33) 
Following the procedure of Sec. Ill, we find 



\Pghz\\v = llPilIWi = 2 ' 



while for the corresponding product operator. 



I I I 

IPghzW'd = ■ 



From here, 

s{pGHz) = {N-l)log2. (34) 

Hence, for the logarithm to the base 2, one gets e{pGHz) = N—1. As it should be, there is no entanglement production 
for a single particle, when N = 1, and the case of two particles, when N = 2, reduces to measures (28) or (31). 



D. Multicat states 



The multiparticle Schrodinger cat states, or, for short, multicat states, are sometimes also called the generalized 
GHZ states. They have the form 

|MC > = ci|ll...l > -h C2I22...2 > , (35) 

where N two-dimensional parts are assumed and the coefficients are arbitrary complex numbers satisfying the nor- 
malization |cip -|- |c2|^ = 1. The related density operator is 

PMC = \MC >< MC\ {TthPmc = 1) . (36) 
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For the norm. (13), we find 

WPMcWv = WpWWi = SUp{|ci|^, |C2|^} . 

The entanglement-production measure (20) becomes 

e{pMc) = (l-7V)logsup{|ci|2, |C2|2} . (37) 

Its value lies in the interval 

0<£(/5mc)< (^-l)log2. (38) 

The maximal entanglement production occurs when |cip = |c2p = 1/2, which goes back to the GHZ states. And 
the entanglement production disappears if |ci| = 1 for any i = 1,2. 

Multicat states can be realized for systems of particles with two internal single-particle states, such as trapped ions 
subject to the action of resonant laser beams [45,46] or Bosc-condensed neutral atoms with an effective interaction due 
to coherent Raman scattering [47]. Instead of internal single-particle states, one can create collective nonlinear states 
by invoking the resonant excitation of topological coherent modes in trapped Bose-Einstein condensates [48-50] . Such 
two-level or two-mode states are usually accompanied by atomic squeezing [50,51]. 



E. Multimode states 



A natural generalization of the multicat states are the multimode states describing a system of N parts, each of 
which can be in one of m different modes. Such a state reads 

\MM > = ^ Ca\nn...n> , (39) 



where 



m 



The corresponding density operator is 

PMM = \MM >< MM\ {TvnPMM = 1) • (40) 
For this case, the entanglement-production measure writes 

£(/3mm) = (1 - A'')logsup]c„]^ . (41) 

n 

It varies in the range 

0<eipMM)<iN-l)logm.. (42) 

The maximal entanglement production is reached when ]c„] = 1/m, while entanglement is absent if any of ]c„] = 1. 
For the two-mode case, when m = 2, one returns to the multicat states. In general, states (39) are related to coherent 
states with a fixed number of particles [40,52,53]. 



F. Hartree-Fock states 



These states, typical of indistinguishable particles, have the structure 



> = 4Fr E |12---^>, (43) 



\HF 

sym 



where a symmetrized or antisymmetrized sum is assumed according to either bosons or fermions are considered. Note 
that here it is implied that A'' particles are in N different single-particle states. The density operator is 
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PHF = \HF >< HF\ {TiH PHF = 1) . (44) 

Calculating 

||pi/F||x> = , \\Pl\\Hi=J^, 

we find the entanglement-production measure 

e{pHF) = log — . (45) 

For N = 2, this reduces to 

s{phf) =log2 (iV = 2), 
as it should be as far as state (43) reduces to the EPR state (26). And for large A'', one gets 

£(pHir) ~ log e (A^^oo). 
It is worth noting that in all examples considered above the limit 

lim £{pn) < oo 

exists. The existence of such a limit is sometimes required as a prercqTiisitc property of any entanglement measure. 
In our case, the existence of this limit is guaranteed by the property of additivity (24). 



V. MIXED STATE ENTANGLEMENT 

Mixed states of physical systems are characterized by density operators that are not idempotent, so that ^ py;. 
Since the entanglement-production measure (20) is defined for arbitrary operators, there is no principal problem of 
applying this definition to any density operators, including the mixed-state density operators. When a statistical 
operator is entangling, the related measure of produced entanglement must be nonzero. 

To concretize the aforesaid, let us consider the operator 

p=i(|ll><ll| + |22><22|) . 

This form corresponds to what one terms a separable statistical operator. However, despite its simple form, the 
operator p produces entanglement, i.e., it is entangling. This is easy to demonstrate by taking, for instance, a product 
function 

/ = 72 ( |1 > +|2 > ) (|1 > +|2 > ) 
from the disentangled set T). The action of p on this function gives 

/5/=^ (|11>+|22>), 

which is a Bell state, that is, a maximally entangled state. Here p really does produce entanglement. Therefore, the 
related measure of produced entanglement (20) has to be nonzero. The latter can be easily verified following the 
general scheme. Thus, we have 

Pi = Trw,^, p = ^ ( |1 X 1| + |2 X 2| ) , 
so that p® = p\® Pi- Calculating the norms 

\\p\\t> = \. \\Pi\\-H, = \, \\P®\\V = \. 

we find £(/3) = log 2, as it should be for an entangling operator p. 
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A. Quantum mechanics 



Let us consider mixed-state density operators obtained by taking partial traces of pure-state quantum density 
operators. For instance, let ns consider the Hartree-Fock pure state density operator (44). Its partial traces produce 
the mixed-state reduced density operators 

P^HF = Trwj,+i • • • T^HnPhf , (46) 
with p = 1,2, . . . N — 1. It is easy to find the norm 

11^ II,,- (^-P)' 
from where the entanglement-production measure (20) becomes 

.(^„,)=,o,(^^^ ,47) 

For a large number of particles and finite p this gives 

^p'hf) ^ ^^^^ loge (TV ^ oo) . (48) 

Measure (47) describes the amount of entanglement produced by the reduced density operator (46) defined for p 
particles from a quantum ensemble of A'' particles. 



B. Statistical mechanics 



Statistical properties of multiparticle systems are characterized by reduced density matrices [41]. Therefore, it 
is natural to consider entanglement realized by these matrices. The general scheme of measuring entanglement 
production can be as follows. 

Let X = {x} be a characteristic space of physical coordinates x, whose concrete nature can be arbitrary. For 
example, x can be a set of spatial Cartesian coordinates, or it may be a set of momentum variables, or a set of 
quantum numbers or mode indices. Assume that there exist several such characteristic spaces Xi = {xi}, with 
i = 1,2, .. .p. And let each space Afj be measurable, with a differential measure dxi allowing for defining the Lesbegue 
integration over X^. In the case of discrete variables, the differential measure is atomic, so that integration reduces to 
summation. The total characteristic space is the direct product 

X=xP^,Xi {p=l,2,...,N) . 

The set x^ = {xi,X2, ■ ■ ■ ,Xp} is an element of X = {x^}. The space X is measurable, with a differential measure 
dxP = lltidXi- 

Let the elements of the single-partite Hilbert space Hi be the vectors ipi = [(fii{xi)] treated as columns with respect 
to the variable Xi, with the scalar product 

vtv'i = i'Pu'P'i) = j V>*i{xi) V>'i{xi) dxi . 
A p-order reduced density matrix is defined as a matrix 

Pp = [Pp{xP,yn] (49) 

with respect to x^ and y^, whose elements are 

Pp{xP, yP) = Tr^V(a;i) • • • V(a;p)pV^(2/p) • • • V'^(2/i) , (50) 

where ipix) is a field operator, p is a statistical operator, and the trace is over the Fock space J^. Definition (50) can 
equivalently be written as the statistical average 

Pp{xP,yP) = <^^{yp)...^\yi)^{xi)...i){xp) > . 
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By this definition, the matrix pp is self-adjoint semipositive. An example is a first-order density matrix 

Pi = [Pi{xi, Vi)] {xi,yi e Xi) , 

with 

pi{x,y) = < tp\y)tjj{x) > . 

A common name for p\ is the single-particle density matrix, which, as was emphasized above, does not imply a 
concrete physical particle but rather tells what is the order of the reduced matrix. 
The trace operations for the density matrices are given by the expressions 



and similarly 



which yields the relation 



T^-HiPi = ^ < ni\p\\ni > = pi{xi,Xi) dxi = N , 

Trwpp = j pp{x^, x^) dx^ = (^^^ , (51) 



. {N-p)\ 



(52) 



For the product operator (19), we now have 



whose norm (14) reads 



To effectively calculate the norms, entering the entanglement-production measure (20), we need to specify the 
single-partite spaces Hi. Recall that the latter can, generally, be chosen in an arbitrary way. However, for density 
matrices, there exists a natural choice related to the eigenvectors of p\ . These eigenvectors are termed natural orbitals 
[41]. Therefore, under jrij >, we shall imply the eigenvectors of p\. And a closed linear envelope of these natural 
orbitals gives the natural single-partite space Hi. Then one can write 

P\ = ^n,n, l^i >< nil ■ 

Keeping in mind the spectral norm, one has 



\n, = sup|L»i.„J . 
A p-order density matrix can be presented as an expansion 

Pp= Y. ^{m,n,} \mi...mp><ni...np\ , (54) 

{rmni} 

SO that 

WPpWv = sup |I?^„^„^}| . 
{"*} 

In the case of a system of identical particles, when all Hi are just copies of the same Hi , then 

llPilIwi = IIpiIIwi • 

Therefore, the entanglement-production measure (20) becomes 

. . _ iN-py.NP\\pp\\^ 

^^PP'-^""^ Nl\\pi\\P,^ ■ ^^^> 

Thus, the problem of quantifying the entanglement generated among any p particles from a given ensemble of N 
particles is reduced to calculating the norms of reduced density matrices. 
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C. Evolutional entanglement 



The reduced density matrices (49), in general, depend on time, 

Pp(^) = bp(a^^y^^)], (56) 

which enters definition (50) either through the field operators i^{x,t) or through the statistical operator p{t). Con- 
sequently, the entanglement-production measure (55) is, generally, also a function of time. Temporal dependence of 
entanglement is named evolutional entanglement [4]. 

As an example of entanglement generated in a nonequilibrium system, let us consider the case when the reduced 
density matrices have the structure of the mixed multimode state 

N\ ^ 

" {N -p)\ Z^^"*^*^ \n...n><n...n\ , (57) 
where Wn{t) are the fractional mode populations with the properties 

0<U'„(i)<l, Y^Wn{t) = l. 

n 

Such multimode, or in the simplest case two-mode, states can be created in trapped Bose-Einstein condensates 
(see reviews [54-56]). This can be done, e.g., by separating a Bose condensate in a two- well, or a multiwell potential 
[55,56]. Another possibility is by generating topological coherent modes in a trapped condensate, by means of resonant 
alternating fields [48 50]. The mixed atom- molecule condensate can also be treated as a two- mode coherent system 
[57]. Coherent collisions between matter waves result in the formation of an effective multimode system [58]. A similar 
system can also be created in the process of superradiant scattering of an atomic Bose-condensed cloud [59-61]. A 
multimode Bose-Einstein condensate has many analogies with coherent optical systems and, in particular, with lasers 
[50,56,62]. 

For the density matrix (57), we get the norm 

"^^^^^"^ " {N-p)\ ' (^^^ 

while for the first-order matrix (52), we have 

M{t)\\n, =N supw„(t) . 

n 

Thus, for the product operator (53), we find 

l|P?(*)l|p= (jyj^^, sup wl{t). (59) 
As a result, the entanglement-production measure is 

£(Pp) = (1 - P) log sup Wn{t) . (60) 

n 

The temporal behaviour of the fractional mode populations Wn{t) are defined by the evolution equations describing the 
corresponding process. For instance, in the case of Bose-Einstein condensates, the evolution equations for the mode 
populations follow from the time-dependent Gross-Pitaevskii equation [48 50,62]. The fractional mode populations 
Wn{t), depending on the physical situation considered, vary in time between and 1. The maximal entanglement 
production happens at that time when all populations coincide, so that Wn{t) = l/m, where m = 1. That is, the 
entanglement-production measure (60) can vary in the interval 

< e{pp) <{p- l)logm . 

Thus, the entanglement-production measure (60) depends on time through the mode populations. The evolution 
of the latter in some cases can be regulated. For example, when the condensate mode structure is due to the 
resonant generation of topological coherent modes, the temporal behaviour of mode populations is governed by the 
applied alternating fields [48-50]. In this way, the time evolution of the entanglement-production measure (60) can 
be controlled, which opens wide possibility for information processing. 
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D. Spin entanglement 



Entanglement production, being defined for arbitrary operators, can be considered for systems of any nature. An 
important class of systems is that corresponding to spin ensembles. Spin entanglement can be studied in the same 
way as entanglement of particles, by invoking a kind of density matrices. 

Let spin operators Sj = {S"}, with a = x,y,z, be associated with a lattice Zjv = {a,} whose lattice sites are 
enumerated by an index i = 1,2, . . . ,N. Similarly to the reduced density matrices (49), composed of field operators, 
we may introduce [40] spin density matrices 



Rp 



d{«/3} 



(61) 



with the elements 



^^<s?:---sf^s'^...s:;> (62) 



composed of statistical averages of spin operators. Form (61) is treated as a matrix with respect to all indices {ij} 

as well as {a/3}, although in principle, one could also consider the matrices Rp"^^ and Rp{ij} defined as matrices 
with respect to only the indices {ij} or {a/3}, under fixed other indices. The first-order spin density matrix (61) is 
i?i = [Rff], whose elements are the correlation functions 

R?f ^<Sfsr>. (63) 
Matrices (61) are self-adjoint, since such are the spin operators. With the traces 

Trn,Ri = Y,Yl < ^1 > = ^^(^ + 1) ' '^nRp = {NSf {S + 1)p , 

i a 

where S is the maximal quantum number of each spin, the product operator (19) becomes 

Rf =Ri®Ri0...^R-^. (64) 

To calculate the operator norms, the single-partite space Tli can be specified as the closed linear envelope of natural 
spin orbitals, which are the eigenfunctions of This is in analogy with the case of reduced density matrices pp. 
Finally, the spin entanglement measure (20) takes the form 

eiRp)=log^^. (65) 
\\Rp \ \t> 

Explicit calculation of this measure will be illustrated in the following sections. 

VI. STATISTICAL THERMAL ENTANGLEMENT 

Instead of considering entanglement realized by reduced density matrices, as in the previous section, one may study 
entanglement produced by a statistical operator p. The latter for a statistical system in thermal equilibrium reads 

p = \ e-P" , Z = Tr e-P" , (66) 

where /? = {kBT)~^, with ks being the Boltzmann constant and T, temperature. One can define, in complete analogy 
with the previous sections, the reduced operators by tracing out some of the states of the total space H. For example, 

pi = T^iH,^.} P ■ (67) 

Entanglement, related to statistical thermal operators, is usually characterized by the concurrence between a pair of 
qubits [63-68] . Here we show that the entanglement-production measure (20) provides a natural characteristic of the 
thermal entanglement, which is produced by the statistical operator (66). 
Let us consider the Ising model 
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N N 

H = --J2 j^j si S^~bJ2S!^ (68) 

ijij i=l 

with an exchange interaction Jij = Jji and an external magnetic field B > 0. Positive interaction Jij > is 
called ferromagnetic, while negative Jij < 0, antiferromagnetic. It is natural to define the basis states \ni > as the 
eigenvectors of 5f . Then the single-partite space Hi is the span of the basis {\ni >}. For concreteness, let us take 
a two-site case, denoting J = J12, and let us consider spin-one-half operators Sf, with S = 1/2. Introduce the 
dimensionless coupling g and magnetic field b by the expressions 



Then we have 



where 



For the reduced operator 



/? J b = (3B . (69) 



P2 = I exp {4<?5f + b {S! + S^)} , (70) 



Z = 2 (e^coshb + . 
p\ = Tr« . p2 , (71) 



we find 



p{ = 1 exp |2ff5f + 6 + I + i exp |-255f + 6 - 



(72) 



The disentangled set T> consists of the states | tT>) I Ti>) I iT> and | ii>. For the norms of p2 and p\, we get 

1 
Z 



P2|h = ^ sup{e«+^ e-»} (73) 



and, respectively, 



\p\\\v = ^ (e^+' + e-^) . (74) 



This, taking into account that pf = pl(^ pi, results in the entanglement-production measure 



S{P2) = log 



2(l + e^^cosh6) , 



(75) 



describing the pairwise spin entanglement production. 

Analyzing the properties of measure (75), it is interesting to compare its behaviour with that of the average 
magnetization per spin 

for which one has 

e^^sinhfo 

" 2(l + e29cosh6) ■ ^ ' 

Expressions (75) and (76) are functions of two variables, coupling g G (—00, -t-00) and magnetic field b G [0, 00). 
In the case of zero magnetic field, when S — > and 6 — > 0, one has M — > and 

elsl 

lini e{p2) = log — — . (77) 
6->o coshgi 
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For low temperature, when T 0, hence /? — > oo and g — > ±00, we get 

lim lim e{p2) = log2 , (78) 

g— >±oo 6^0 

that is, the maximal pairwise entanglement production. In the opposite case of high temperatures, when T ^ 00, so 
that /3 — > and 5 — > 0, we find 

lim lim e(p2) = , (79) 
s— >o 6^0 

which means that high-tcmpcrature fluctuations destroy entanglement production. Note that the limits 6 — > and 
g commute with each other. 

If the external magnetic field increases, with B — > 00 and 6 — > 00, then 

lim e{p2) = , lim M = i . (80) 

This shows that there is no entanglement production between perfectly aligned spins. 

When temperature diminishes, such that T ^ and /? 00, but the magnetic field is nonzero, B ^ 0, then the 
resulting expressions depend on the relation between 6 — > 00 and \g\ — > 00. Thus, we obtain 

log2 {b + 2g -00) 
lim e{p2) = { Zoff(3/4) {b + 2g^ 0) (81) 
(6 + 2g-*+oo). 



And for magnetization (76), we find 



{b + 2g^ -00) 
lim M = { 1/6 {b + 2g^0) (82) 
1/2 {b + 2g^ +00). 



Magnetization plays the role of an order parameter. The above analysis shows that the larger is the magnetization, 
the smaller is the entanglement-production measure. In this case, entanglement production and order are complimen- 
tary to each other. To illustrate this by another example, let us consider the Hamiltonian (68) with ferromagnetic 
interactions Jij > of long-range type, when Jij depends on N so that 



1 ^ 

lim Jij = , lim — > J^, < 00 . (83) 

JV->oo N^oo N ^ 

In that case, as is known [69], the Hamiltonian (68) is asymptotically, as A'' — > 00, equivalent to the mean-field form 
where 



Hi 



<5|> +B^S! + l ^-^'^ <s!><s^> . 



This implies that the statistical operator (66) is asymptotically, as N 00, equivalent to 



But then pN ^ p% and we come to the limit 



lim e{pN) = . (85) 

iV— >(x) 



Long-range ferromagnetic interactions, satisfying condition (83), organize long-range magnetic order in the spin system 
with a finite magnetization M. At the same time, this yields the absence of entanglement production between 
ferromagnetically aligned spins. 
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VII. ENTANGLEMENT AND PHASE TRANSITIONS 

The analysis of the previous section hints that there should be a relation between the entanglement-production 

measure and an order parameter. The latter, in turn, experiences dramatic changes at phase transitions. Hence, 
entanglement production may also exhibit essential changes under phase transformations, which is illustrated below. 



A. Bose-Einstein condensation 



Let us study the entanglement realized by the reduced density matrices (49). The properties of the latter are 
described in detail in book [41]. At high temperature, much larger than the condensation temperature T^, one has 

llPplk^llPill?,, (T»Te). (86) 
Therefore, the entanglement-production measure (55) is 

e{p,) log ^ (T » T,) , (87) 

which is typical of the Hartree-Fock form. 
At low temperature T <C Tc, we find [41] that 

WppW^ ^ -(j;rhy_ ' WpiWh^^n. (88) 

Therefore, measure (55) becomes 

eifp) ~ (T « T,) . (89) 
This means that entanglement production diminishes when the Bose-Einstein condensation occurs. 



B. Superconducting transition 

At temperatures much higher that the critical temperature Tc, density matrices are of the Hartree-Fock type, which, 
for large N ^ 1, yields 

loge (T»T,). (90) 

At temperatures below T^, the structure of the reduced density matrices essentially changes, as is thoroughly described 
in Ref. [41]. Then one has 

II II f Ar(3'-i)/2 ^ ^^^^ 

m4 



where Cp is a constant of order one. From here 



Thus, for T <Tc, finite p, and large N ^ 1, we obtain 

i'^'^'K (92) 
I 2 log ^ (P even) . ^ ' 

Comparing Eqs. (90) and (92), we see that the entanglement-production measure £{pp) increases under arising 
superconductivity. This is contrary to what happens under Bose-Einstein condensation. Such a difference should 
not be surprising and it can be easily understood as follows. Considering here entanglement production, the single- 
partite space Wi has been treated as the space of single-particle quantum states. In the case of the Bose-Einstein 
c;ondcnsation, there appears ordering of particles, which leads to the decrease of their entanglement production. But 
under superconducting transition, the arising order has to do with pairs of particles, that is. Cooper pairs, and not 
with separate particles. If we define the single-partite space as the space of the Cooper pairs quantum states, then the 
overall situation would become similar to that occurring at Bosc-Einstcin condensation. Then the order appearing 
between Cooper pairs under superconducting transition would result in the diminishing entanglement production of 
these pairs. 
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C. Magnetic transition 



To study the interplay between magnetic ordering and entanglement production, we shall invoke the spin density 
matrices (61). As a particular case, let us consider these spin density matrices composed of the z-components Sf of 
spin operators, with i = 1,2, ... ,N enumerating lattice sites. Thus, the spin density matrix Rp = [R^ijy] possesses 
the elements 

Rw}^<Sl---SlS!^---Sf^> . (93) 

Clearly, this is a Hermitian semipositive matrix. 
The eigenfunctions of Ri = [< Sj Sf >] are the vectors 

\k> = [vPfe(ai)] , ^fc(a) ^ ^ e*- , (94) 

V J\ 



treated as A^-order columns with respect to the lattice vectors aj G Zat. Functions (94) form a complete orthonormal 
basis {\k >}. A single-partite space is defined as the closed linear envelope Hi = jC{\k >}. The norm of Rp over V 
can be calculated as 

||-Rp||x> = sup < ki . . .kp\Rp\ki . . .kp > . (95) 
{fei} 

The elements of Ri are the correlation function < > having the property [70] 

lim < > = M^ , (96) 

an — *oo 



where Uij = la, — aj| and M is the magnetization 



1 ^ 



Here wc keep in mind ferromagnetic phase transition. Calculating norms (95), we use, for simplicity, the mean-field 
approximation and assume large N oo. Above the critical temperature Tc, where M = 0, we find 

\\Rp\\T, = {2p-l)\\S^P {T>T,), (97) 

with 

(2'dV 

{2p - 1)!! = 1 X 3 X 5 X . . . X (2p - 1) = . 
And below T^, where M ^ 0, we have 

||i?p||p~Aff M^f (T<Te). (98) 

In particular, 

\\Ri\\h^=S^ + NM^ 

at all temperatures. Therefore, for the entanglement-production measure (65), we obtain 

e{Rp)=log^^_ (T>Te) (99) 

above the transition temperature. For instance, 

£(i?2)=log3, e(E3) = log 15 . 
And below where there appears magnetic order, wc get 

e{Rp) ~ (T < Te) (100) 

for any finite p. 

As is seen, the arising magnetic order also leads to the decrease of entanglement production, similarly to what hap- 
pens at Bosc-Einstcin condensation. The same qualitative relation between entanglement production and ordering for 
different phase transitions can be understood if to keep in mind that a ferromagnetic phase transition is accompanied 
by the condensation of magnons [70]. 
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D. Order indices 



The intimate relation between entanglement production and ordering, illustrated above by several examples, can 

be explained by the most general arguments. For this purpose, wc need to resort to the notion of operator order 
indices that can be introduced for arbitrary operators [40], and, in particular, for reduced density matrices [41]. The 
operator order index of an operator A is defined [40] as 



logjTr A\ 



In describing the ordering in physical systems, the role of operators A is played by reduced density matrices. The 

latter can be the standard density matrices constructed of field operators [41], as well as spin density matrices, lattice 
density matrices, or matrices composed of other operators [40]. The larger is the order index uj{A), the higher is the 
level of ordering corresponding to the operator A. Such order indices characterize both long-range as well as various 
types of mid-range order. They are suitable for describing off-diagonal and diagonal orders and can be applied for 
any physical system, equilibrium or nonequilibrium, infinite or finite. Also, they do not involve the notion of broken 
symmetry and can be employed wlic;n the order parameters are not defined. 

If the operator A represents a density matrix, it is semipositive. Then \ \A\\ < TrA, because of which u!{A) < 1. 
Assume that a semipositive operator A is associated with a kind of ordering in a physical system. The order is absent 
when ||A|| < TiA, and then uj{A) < 1. 

There are two types of long-range order that may develop in physical systems, total and even orders [40,41]. Under 
the arising total long-range order, the norms of the related density matrices increase, so that \\pp\\ ^ Trpp, hence 
(-Op{p) 1. The order is total, which implies that the property ||jOp|| ^ Trpp is valid for the density matrices of all 
orders p = 1, 2, . . ., including ||pi|| ~ Trpi. From here, \\pf\\ ~ Trp®. According to the normahzation condition (18), 
we have Tr pp = Trpf. Therefore, jjppjj ^ [jp® ][, which results in e{pp) — > 0. As is shown above, this situation takes 
place at Bose-Einstein condensation and ferromagnetic transition. Thus, under the appearing total long-range order, 
the order indices increase but the entanglement-production measure decreases, 

io{pp) ^ 1 , e{pp) ^ {total order) . 

Here increasing total order implies diminishing entanglement production. 

The situation is different for even long-range order. Then the norms of the density matrices also increase, but so 
that [40,41] 

^ r ^Tipp/N {p odd) 
\ ^Trpp {p even) . 

For instance, ||pi|| ~ const but ||p2|| ~ iV. Because of this, the order indices are different for p odd and p even, 

'^yPp)^\l/2 {peven). 

There is no order in the single-particle matrix pi, as far as co{pi) = 0, but u>{p2) = 1/2. For all p = 1, 2, . . ., we have 
<jJ2p-i < W2p, and they coincide only in the limit p ^ oo. Such a behaviour is typical of superconducting transition. 
Entanglement production is also different for odd and even numbers of particles, 



^ log TV [podd) 
I log N (p even) 



The entanglement-production measure increases, since the inequality ||pi|| <C \\pp\\, with p > 1, yields || <C \\pp\\- 
Now both the order indices for p > 1 and entanglement-production measure increase together. Hence, increasing even 
order results as well in increasing entanglement production. 

In conclusion, a general definition of entanglement-production measure for arbitrary operators is introduced. The 
concept is valid for systems of any nature. The interplay between entanglement production and phase transitions is 
elucidated. 
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